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The recent Planck results and future releases on the horizon present a key opportunity to address 
a fundamental question in inflationary cosmology of whether primordial density perturbations have 
a quantum or thermal origin, i.e. whether particle production may have significant effects during 
inflation. Warm inflation provides a natural arena to address this issue, with interactions between the 
scalar inflaton and other degrees of freedom leading to dissipative entropy production and associated 
thermal fluctuations. In this context, we present relations between CMB observables that can be 
directly tested against observational data. In particular, we show that the presence of a thermal bath 
warmer than the Hubble scale during inflation decreases the tensor-to-scalar ratio with respect to the 
conventional prediction in supercooled inflation, yielding r < 8|n t |, where nt is the tensor spectral 
index. Focusing on supersymmetric models at low temperatures, we determine consistency relations 
■ between the observables characterizing the spectrum of adiabatic scalar and tensor modes, both 

for generic potentials and particular canonical examples, and which we compare with the WMAP 
and Planck results. Finally, we include the possibility of producing the observed baryon asymmetry 
during inflation through dissipative effects, thereby generating baryon isocurvature modes that can 
' be easily accommodated by the Planck data. 
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I. INTRODUCTION 



The recently announced results of the Planck satellite [26| seem to paint a simple picture of the early universe, with 
the observed temperature anisotropies in the Cosmic Microwave Background (CMB) being, at least on small scales, 
well described by a primordial spectrum of nearly scale invariant, gaussian and adiabatic fluctuations evolving in a 
ACDM cosmological background. Although these can be generated by a single slowly-rolling scalar field dominating 
the early stages of an otherwise empty universe, one may wonder whether this really satisfies Occam's razor, in 
I/"") ' particular taking into account that the simplest monomial potentials in chaotic inflation are now either ruled out or 
are barely allowed by Planck. Is a single scalar field rolling down a nearly flat potential, which is either fine-tuned or 
' protected by a special symmetry, and with little interactions with any other degrees of freedom the simplest description 
, of the observational results? 

A fundamental question to address is whether the data really establishes an inflationary universe in a true vacuum 
state or if more general statistical distributions may yield equally good or even better descriptions of the observational 
results. A particularly simple and natural possibility is the presence of a thermal state during inflation, which despite 
the accelerated expansion rate may be sustained if for example the inflaton itself can act as source of particles that 
quickly thermalize. Although the details of particle production may differ depending on the nature and structure of 
the interactions between the inflaton and other degrees of freedom, this generic inflationary paradigm is known as 
warm inflation [H-ll (see also 0-0; in an observational context see [10H12| ). with the dynamics of the background 
field and associated perturbations being governed by fluctuation-dissipation dynamics. 

Warm inflation thus provides a natural background to test alternatives to the standard inflationary paradigm 
and address the general question of whether primordial density perturbations have a quantum or thermal origin. 
The announcement of the Planck results and the prospects for new data releases concerning CMB polarization thus 
motivate a closer look at this issue, in particular determining which observables and relations between them are 
suitable to address this question in the light of the new and future observational results. 

Warm inflation can be successfully realized through a generic supersymmetric two-stage model where fields that 
couple directly to the inflaton acquire a large tree level mass due to the inflaton's large vacuum expectation value 
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(vev), but may subsequently decay into light degrees of freedom, which upon thermalization create a thermal bath 
concurrent with the inflationary expansion (13|. If the inflaton vev is sufficiently large, the masses of these mediator 
fields become heavier than the temperature of the thermal bath, which is known as the low-temperature regime, 
and thermal corrections to the inflaton potential are thus Boltzmann-suppressed, preserving the tree-level flatness of 
the potential. Despite supersymmetry being broken due to the finite energy density and temperature during warm 
inflation, the supersymmctric nature of the model keeps both radiative and thermal corrections to the potential under 
control. This mechanism leads to dissipative effects that further damp the scalar inflaton's motion, allowing for longer 
periods of accelerated expansion, and the continuous production of relativistic particles may automatically yield a 
graceful exit into the standard radiation-dominated cosmology. 

As proposed in [T3|) if the heavy mediators decay through B- and CP- violating interactions, the production rates 
of baryons and anti-baryons will be different and a baryon asymmetry may be produced during inflation, a mechanism 
known as warm baryogenesis. Dissipation is naturally an out-of-equilibrium process, with annihilations of the light 
degrees of freedom being inefficient in repopulating the classical background condensate, thus satisfying the Sakharov 
conditions. The structure of the interactions closely resembles that of GUT baryogenesis models, similarly requir- 
ing at least two heavy mediators with different masses and complex Yukawa couplings, and in fact the two-stage 
superpotential is a natural feature of GUT- like gauge theories and associated D-brane constructions [l5j |. 

The asymmetry approaches a steady state solution during inflation, due to the slowly varying temperature of 
the thermal bath and will freeze out once inflation ends and the Hubble parameter drops sufficiently to make the 
dissipative mechanism inefficient. The produced asymmetry is suppressed in the low-temperature regime, allowing for 
0(1) couplings with CP violating phases and mass differences at the few percent level whilst producing an asymmetry 
within the observed window. This is in contrast to standard thermal GUT scenarios where heavily suppressed couplings 
are required to generate the observed asymmetry. Furthermore, standard thermal GUT baryogenesis models require a 
reheating temperature above the GUT threshold in order to generate a baryon asymmetry from an initial population 
of heavy fields in thermal equilibrium, which may lead to an overproduction of dangerous GUT relics or gravitinos. 
On the other hand, the temperature during warm inflation must be below the heavy mass threshold, preventing the 
production of GUT monopoles and potentially solving the gravitino problem for sufficiently strong dissipation [16[ . 



As the baryon asymmetry produced through dissipation depends on the inflaton vev and the temperature of the 
thermal bath, the resulting baryon-to-entropy ratio will exhibit fluctuations that are fully correlated or anti-correlated 
with the adiabatic inflaton fluctuations and induce baryon isocurvature perturbations that are subsequently imprinted 
on the CMB. Warm baryogenesis thus avoids the problems that plague standard models of thermal GUT baryogenesis 
whilst simultaneously being a testable model of baryogenesis. 

In warm inflation, the observables that characterize the density and tensor perturbations spectrum depend not only 
on the slow-roll parameters, as in cold inflation, but also on the value of the dissipative coefficient when the relevant 
scales exit the horizon during inflation, Q = T/3H, and on the number of relativistic degrees of freedom, g*. In 
general, this yields too many parameters to allow one to fully relate the observable quantities. On the other hand, 
the generation of baryon isocurvature perturbations provides an additional observable and, furthermore, canonical 
potentials used in inflationary model building are characterized by only two or three parameters. As we will show, 
this will allow us to derive several consistency relations that can be directly tested against observational results. 

In this work, we focus on consistency relations between observables based on the above mentioned two-stage 
supersymmetric warm inflation model, in the low-temperature regime, assuming the observed baryon asymmetry is 
produced during inflation and that adiabatic perturbations originate solely from thermal fluctuations of the inflaton 
field. Here, we anticipate the main results of this work, which are the consistency relations for weak and strong 
dissipation, Q -C 1 and Q > 1, respectively, that hold for any generic inflationary potential. These relate the scalar 
and tensor spectral indices, n s and rij, the associated running parameters, n' s and n' t as well as the running-of-the- 
running of the scalar index n", the tensor-to-scalar ratio, r, the fraction of baryon isocurvature perturbations, Bb 
and its spectral index, rii SO : 
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We will place these relations in context below and apply them to some canonical potentials. These relations correspond 
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to the dissipation coefficient obtained in the low-temperature regime through the excitation of heavy virtual modes 
coupled directly to the inflaton, which is the most well-studied case in the literature (see e.g. 0, [lg, [TtJ for recent 
studies of warm inflation dynamics in this scenario). In this sense, these results cannot be used to probe the full warm 
inflation paradigm. One should note that, in the low-temperature regime, the excitation of on-shell models may give 
the dominant contribution to the dissipation coefficient for sufficiently small Yukawa couplings (l4| . a scenario that 
is, however, outside the scope of this work. 

On the other hand, as we show below, in warm inflation models we find: 

r < 8|n t | , (2) 

independently of the form of the dissipation coefficient, with r = 8|n t | corresponding to cold single-field models. This 
is also a generic feature of multiple-field inflationary models [l8[ , so that the consistency relations in Eq. (JIJ may be 
crucial in distinguishing these from warm inflation models if r < 8\n t \ is observed. 

We begin by introducing the relevant dynamics of warm inflation in Section [Til before going on to look at the 
relevant observables in Section IIIII We then analyze in more detail a few canonical examples, namely monomial, 
hybrid and hilltop models, in Section llVI In Section [V] we discuss our main results and the prospects for the future 
releases of CMB data by the Planck collaboration. 



II. WARM INFLATION DYNAMICS 



Warm inflation can be realized in supersymmetric theories through a "two-stage" generic superpotential of the 
form: 

W = |$x 2 + *±XY? + /($) , (3) 

where the scalar component of $ yields the inflaton field and /($) is a general holomorphic function describing the 
self-interactions in this sector. The X components couple directly to the inflaton and become heavy, while the fields 
in the Y sector remain light and form the radiation bath. The coupling of the inflaton to heavy fields which decay 
into light degrees of freedom then induces dissipative particle production, yielding an additional source of friction in 
the inflaton equation of motion and a source of quasi-thermal radiation: 

4> + (3H + T)0 + = , p R + AHp R = , (4) 

where indicates the derivative of the potential with respect to the inflaton field, T is the dissipative coefficient 
and pn the radiation energy density. It has been recently pointed out that bulk viscous effects in the quasi-thermal 
bath may modify the dynamics and prolong inflation [19[ , although we will not consider this case in our discussion. 
Accelerated expansion occurs in the slow-roll regime, where the equations of motion reduce to: 

3H(1 + Q)4> » -V* , 4 PR w 3Q0 2 , (5) 

with Q — T/3H as defined above. In warm inflation, slow-roll is achieved if the following conditions are satisfied: 

e * = ! f (y) 2<1 + g ' ^ = ™p(it) <! + Q> = < l + Q> ( 6 ) 

where m p denotes the reduced Planck mass. Other slow-roll parameters ensure that the radiation is not redshifted, 
that dissipation does not cause the radiation to dominate too soon and that thermal and quantum corrections to the 
potential are not too large (see e.g. Q). In the low-temperature regime, the form of the dissipative coefficient depends 
on whether the mediator fields in the X sector are produced dominantly on- or off-shell [14j . We will focus on the 
latter case, although the production of on-shell modes gives a promising avenue for model building that has yet to 
be analyzed in detail. The dominant source of dissipation then corresponds to the production and decay of virtual 
scalar modes [3, H(| , giving: 

T « , (7) 

where the constant C# depends on the field multiplicities and couplings in the light and heavy sectors [l4| • From the 
slow-roll equations we can relate Q and <j) via: 
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Using these relations, the evolution of the dynamical quantities can be written in terms of the slow-roll parameters: 

<f> (i + Q)' [> 

Q'= (1+ Q 7Q) (lO^-6^ + 8a ) , (10) 



fry (i + 7Q)Vi + O v " i + Q 7 U 

where the prime indicates a derivative with respect to the number of e-folds, dN e = Hdt. Note that Eq. (J9)) is 
independent of the form of the dissipative coefficient. 

For T > H, the dominant contribution to the primordial perturbation spectrum are the thermal fluctuations from 
the coupled inflaton-radiation system. Upon exiting the horizon these thermal fluctuations freeze out as classical 
perturbations and during slow-roll the amplitude of the curvature perturbation power spectrum is given by 0, [2l| : 

where all quantities are evaluated at horizon-crossing. Note that the cold inflation limit is obtained when T/H — > 1 
and Q — > 0, i.e. in the absence of dissipation the temperature corresponds to that of the cosmological dc Sitter horizon. 

III. OBSERVABLES 

The dependence of the curvature power spectrum, Eq. (TT2")) . on the scale of the perturbations is given by the scalar 
spectral index 

n s - 1 = ^ = (1 + Q) ; i + ?Q) (-(2 + 9Q)e, - 3Q V , + (2 + 18Q W ) . (13) 
The running of the scalar spectral index is given by: 

+ (1 + QK1 + 7Q)2 (5 ^ " ^ + 4fT0) ' (14) 
where Aq = Q/(l + 7Q) and primes denotes derivatives with respect to the number of e-folds, in particular: 

£ = (I + Q) ( 2 ^ 1 ^ = (1 + Q) ~ ' ^ = (1 ^Q) ^ + 2 ^ ~ ^) ' ( 15 ^ 

with £4, = 2m 2 p {V < j >< j >< j > /V^). 

Tensor modes are not affected by dissipation due to the weak coupling of the graviton to the thermal bath and so 
the tensor perturbation spectrum is governed by quantum fluctuations. The tensor-to-scalar ratio is then modified 
with respect to cold inflation by an additional thermal factor: 

H\ 16e 

We note that as long as the inflaton is solely responsible for the generation of adiabatic perturbations, which are 
thermally produced, we may use the slow-roll equations to show that: 



7T pR. \ m 



H ,2 



(17) 



which is analogous to the cold inflaton result. However, one should note that thermal fluctuations are generically 
larger than quantum fluctuations, yielding the same amplitude for different field values. In particular, in large field 
models, one may obtain the observed amplitude for <ft < m p pj. 
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The scale dependence of the tensor power spectrum is given by the tensor spectral index, which due to the decreasing 
Hubble parameter is always negative: 

d In Ph 2eA, 

Ut = "dlnT = ~T~K? ~ ' (18) 

where nt ^ —2 in order to satisfy the slow-roll conditions. In cold inflation with a single slowly rolling scalar field 
we thus find r = 8|nt|, which is a well-known consistency relation. In warm inflation, the relation between the 
tensor-to-scalar ratio and the tensor spectral index is, on the other hand, given by: 

for T > H, while one recovers the equality for T/H — > 1 and Q — > 0. Thus, if a future experiment measures r and rit 
and finds r < 8|n t |, this could be a hint for a thermal origin of the adiabatic perturbations, although this is also the 
case for generic multiple field models as mentioned earlier and other observables must be used in order to break this 
degeneracy. 

There is also the possibility that the observed baryon asymmetry is produced during warm inflation through 
dissipative effects as described earlier. In this scenario, the fluctuations of the inflaton are directly transferred to 
the baryon-to-entropy ratio which yields baryon isocurvature perturbations that freeze out during inflation and are 
imprinted on the CMB j22j. These isocurvature perturbations, given by Sb — 5rj s /rj S: are thus fully correlated or 
anti-correlated with the adiabatic perturbations generated by the inflaton field. It is convenient to express this in 
terms of the ratio between the baryon isocurvature perturbations and the curvature perturbation which for the 
dissipative coefficient in Eq. J7]) is given by [f4| : 

B b = y= (1 + q) 2 2 (1 + 7Q) ( 2 ^(1 + Q) - e*(3 + Q)~ + 5Q)) . (20) 

Whether isocurvature perturbations are correlated or anti-correlated with adiabatic perturbations depends upon the 
sign of Sb. 

The spectral index of the isocurvature perturbations is in principle another observable that can be used to test the 
warm baryogenesis scenario. For the dissipative coefficient in Eq. ([7]) the isocurvature spectral index is given by: 

dlnS B 7Q 

2e 4> (2r ]4> - ^)(1 + Q) - 2e (2e - ry )(3 + Q) - ct (ct + 2e - ^)(3 + 5Q) 



+ 



(I + Q)(2^(l + Q) - e (3 + Q) - (i (3 + 5Q)) 
(2?70 - €0 - 5<70) (IQe^ - 6?70 + 8a$)Q 2e < j > - ^ 

(1 + 7Q)(27fo(l + Q) - e (3 + Q) - a (3 + 5Q)) ~ (1 + Q) 
e (2 - Q) - 714,(2 + 5g) + 2a (l + Q) 
2(1 + 7Q)(1 + Q) 



(21) 



This implies that the baryon isocurvature spectral index will typically be proportional to a combination of slow-roll 
parameters for weak dissipation, which is further suppressed by the dissipative ratio Q for strong dissipation (see 
below). Hence, we expect this to be small in general, although formally it can take any positive or negative value 
depending on the form of the inflaton potential. For comparison, in the curvaton scenario the spectral index of the 
(cold dark matter) isocurvature spectrum is the same as for the adiabatic spectrum (see e.g. (23[), whereas isocurvature 
modes produced from the quantum fluctuations of axion fields correspond to an almost scale invariant spectrum [2~i| . 

The bounds on anti-correlated cold dark matter isocurvature modes from the 9-year WMAP data analysis (WMAP 
9), using the combined WMAP+eCMB+BAO+f/o data and taking into account that Hcdm/^b ~ 5, yield an upper 
bound \Bb\ < 0.33 for Bb < [23. Fully correlated modes have not been analyzed with WMAP 9 data, with an 
earlier analysis yielding Bb < 0.43 [23j |. 

In its recent data release, Planck has placed new bounds on the presence of a combination of cold dark matter and 
baryon isocurvature modes (CDI), which are indistinguishable from the CMB point of view. For a general admixture 
of adiabatic and isocurvature perturbations, independently of their correlation or the latter's spectral index, the 
Planck results yield /3j so < 0.075 for the CDI fraction at the comoving wavenumber fci ow = 0.002 Mpc" 1 [26|. This 
translates into a bound \Bb \ 1-5, which is surprisingly less constraining than the WMAP results. According to the 
Planck collaboration, this is related to the fact that the data actually prefers models with a significant contribution 



6 



of CDI (or neutrino density isocurvature modes), driven mainly by the amplitude deficit observed at low multipoles. 
Although no definite evidence for isocurvature modes was found in the data, it nevertheless allows for a substantially 
large fraction of isocurvature modes. 

Given the model-dependence of the sign of B B and the value of rii SO discussed above, in our analysis below we 
will include the general bound from Planck and the earlier WMAP bounds to illustrate our results. One should 
nevertheless keep in mind that, for particular potentials, rii SO will be related to n s — 1 and more stringent bounds 
should apply as for the axion and curvaton cases. For example, for monomial potentials, V(4>) oc 4> n 1 and a dissipative 
coefficient of the form Eq. ([7]), we find: 



n s — 1 2 n — 2 ' 
n iso _ 42 - lln 
n s - 1 ~ 84 - 30n ' 



<2<i , 

Q > 1 • (22) 



We would like to emphasize that the warm baryogenesis scenario may generate isocurvature perturbations which are 
either correlated or anti-correlated with adiabatic perturbations, depending on the inflaton potential, whereas the 
recent literature has focused mostly on anti-correlated and uncorrelated modes, motivated by curvaton and axion 
models. We hope that this work motivates the need to look for fully correlated isocurvature modes in future data. 

In general, for the low-temperature dissipative coefficient in Eq. ([7]), the expressions for the observables are too 
complicated to yield simple relations between them. If, however, we look at the strong (Q ^ 1) and weak (Q <C 1) 
dissipative regimes separately, these expressions simplify considerably and such consistency relations can be obtained. 
To illustrate these relations, we will use the 68% and 95% C.L. contours in the (n s , n' s ) plane obtained from the Planck 
data for a ACDM cosmology with a running scalar spectral index with and without tensor perturbations [26| . 

1. Weak dissipation 
In the weak dissipative limit the observables reduce to: 



n s - 1 « -2e + 2a , n' s w -2e' + 2a' , < « -2e'i + 2a'l , r 



16e / 16n 2 P n C R \ 1/3 



1/3 V 3 J 



u ofo o o n 94 + 2e*(&-4ifc+30>)+<7*(60*-5rfe) 

n t « -2e , B B « 2(2^ - 3e - 3ct ) , n iso « — v - — . (23) 

3e^ - 2^ + 3o> 

From these expressions we can derive two consistency relations for an arbitrary inflationary potential with the low- 
temperature dissipative coefficient in the weak dissipative regime, Q <C 1: 

B B =An t -9(n s -l) + 8 ^-^- 4 < , (24) 

n s — 1 + n t 

_ _ 1 ( 5{ns _ 1} _ 4nt , 8 (™s " 1)K - 2(n a - l) 2 ) , 6((n s - l) 2 - 2n' s ) 



''ISO 



4 1 ' (n s -l + n t ) 2 n s -l + n t 

-7 + 7nl + 16^ - 3n t - 12<(3 + n t ) - 3n^(7 + n t ) + 3n s (7 + 12<. + 2n t ) 
l + n 2 s + An' s - n t (5 + 4n t ) + n s (5n t - 2) 



(25) 



Figure Q] shows the contours of B B from Eq. (|24p in the (n s ,n' s ) plane for n t = —0.01. For < 0.01 the effect 
of the tensor index becomes negligible in Eq. (1241) and the parameter space plot tends to this limiting form. For 
n t < —0.4 we do not find any solutions within the 95% C.L., so that this places an upper bound on \n t \ in the weak 
dissipative regime in the warm baryogenesis scenario. 

2. Strong dissipation 



In the strong dissipative regime there is no consistency relation involving only n s ,n' s ,r,n t , B B , as there are 5 
parameters (e^, 77^, o^, Q) for a given value of and only 5 observables. In the strong dissipative limit these are 
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FIG. 1: Contours of Bb in the (n 3 ,n' g ) plane for Q <C 1, with n t = —0.01. The dark grey excluded region represents the current 
bounds on Bb from Planck, while the lighter gray indicates previous limits from WMAP. 



given by: 



n„ ~ 1 » 7Q ( 6cr ^ ~ 3e ^ ~ ' w ~49Q2 ( 54e + 6 ^ ~ 2 ^ a 4> + 6o > ~ ^(7^ + 20^ + 48o- ) , 

r ^Q^i 3 J ' n *"- 2 Q' S fl «— (2^-e,-5^). (26) 

However, if we somewhat forward-thinkingly include the running of the tensor spectral index: 

< = - 2 TTq + 2 (TTW = iTTwhW) { - {2 + 9Q)e * + * + 4Q( ^ + **» ' (27) 

we obtain the following consistency relation: 



1152B 2 B n* {16TT 2 P n C R \ 1/3 



(n t (8(l - n s ) + 3n t ) - 3n' t ) 2 \ 3 / 



(28) 



We note that the dependence of this relation on the number of relativistic degrees of freedom, g*, is rather mild. For 
example, for = 0(10 — 100), the tensor-to-scalar ratio varies only by a factor of ~ 2. 

If in addition, we may consider the isocurvature spectral index, which in the strong dissipative regime is given by: 

_ 27e 2 + 28e^ - 17e ?7 - 18r$ - 3e cr - 29^^ + 80a 2 
niS ° ~ 14Q(e - 2^ + 5a ) ' (29) 

so that we can derive the following consistency relation relating rii so , n s , n' s , n t and n' t : 

_ n\ [-32((n s - l) 2 + 4<) + 4(n s - l)n t + 9nf\ + 2^(2(1 - ra s ) + Ahn t )n' t - 3n: 



8th [8n t (l - n s ) + 3n 2 - 3raJ] 



/2 

(30) 



However, since we do not expect the running of the tensor spectral index nor the baryon isocurvature spectral index 
to be accurately measured in the near future, we will not include the consistency relations in Eq. (|25p . Eq. (|28[1 and 
Eq. ([50)1 in the remainder of our analysis. 



IV. INFLATIONARY MODELS 



In the previous section we have derived general consistency relations between observables that are independent of 
the form of the inflationary potential in both the weak and the strong dissipative regimes. If we consider specific 
potentials, the slow-roll parameters are typically no longer independent and we may find stronger consistency rela- 
tions. In particular, we obtain, in some cases, consistency relations that are independent of the baryon isocurvature 
perturbations, so that these relations hold even if the baryon asymmetry is not produced during warm inflation. In 
this section, we thus consider a few canonical models, in particular monomial, hybrid and hilltop potentials. 
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A. Monomial potentials 



Monomial potentials have the generic form V — Vo(<p/m p ) n . In cold inflation, one typically requires super-planckian 
field values in order to achieve ~ 40 — 60 e-folds of accelerated expansion. For warm inflation, on the other hand, this 
can be achieved in the sub-planckian regime due to the additional friction. 

The slow-roll parameters for the monomial models are given by: 

£ '4(t) ! ' % = »<»-l>0) 2 . ».= <SV f* = ^-l)<»-2)(^)\ (3D 

This yields a blue spectral index, n s > 1, for the quadratic potential (in contrast to cold inflation), which is 
disfavoured by the latest Planck data. On the other hand, for the quartic potential one finds a spectral index in 
excellent agreement with the Planck measurement, even in the weak dissipative regime Q. 



1. Weak dissipative regime 

In the weak dissipative regime, we derive the following consistency relations: 

Bb = (1-».)-^, »,-(«.-l)- 5 Ji 3J . (32) 

Note that the tensor spectral index is always negative, so that these consistency relations only hold in the regions 
of the (n s ,n' s ) plane where n t < 0, i.e. the model can only yield values of the scalar spectral index and running in this 
region. Also, nt > —2 in order to satisfy the slow-roll conditions when the observable scales leave the horizon during 
inflation. 



Baryon Isocurvature fraction Tensor spectral index 




0.92 0.94 0.96 0.98 1.00 0.92 0.94 0.96 0.98 1.00 



FIG. 2: Contours of Bb (left) and nt (right) in the (n s ,n' s ) plane for monomial potentials in the weak dissipative regime. We 
exclude the grey (light grey) regions where Bb is outside the observational bounds for Planck (WMAP) (left) and where the 
consistency relation yields n t > or n t < —2 (right). 

As we illustrate in Figure [2j for n s < 1 the allowed regions where Bb is within the observational bounds and 
—2 < n t < overlap with the 95% C.L. region, with n t falling within the 68% C.L. contour for n t ~ —0.1. Also, 
a positive running would indicate a positive correlation between the isocurvature and adiabatic perturbations with 
Bg ^ 0.1 and a tensor spectral index n t > —0.1. A negative running may, on the other hand, give anti-correlated 
isocurvature modes. 



2. Strong dissipative regime 

In the strong dissipative regime the value of Q at horizon-crossing becomes an additional parameter determining 
the relevant observables compared to the weak dissipative regime. Thus, we may only obtain a single consistency 



9 



relation: 



(33) 




FIG. 3: Contous of n t in the (n a ,n' s ) plane for monomial potentials in the strong dissipative regime. We exclude the grey region 
where the consistency relation yields n t > or nt < —2. 

In the strong dissipative regime, nt is consistent with the 95% C.L. limit and only for larger values, nt ~ —0.5, 
does nt fall within the 68% C.L. contour, as shown in Figured] For n s < 1 (> 1) the running is necessarily negative 
(positive) , which gives a strong constraint independent of whether the observed baryon asymmetry is produced during 
warm inflation. In this regime, there is no consistency relation involving Bb, but since Eq. (J33J) only holds for Q 3> 1 
this effectively limits the value of Bg at each point in the (n s ,n' s ) plane. 



B. Hybrid potentials 

We will now consider an example of small field models, hybrid models, where the inflaton couples to one or more 
'waterfall fields' with a Higgs-like potential. These fields couple to the inflaton in an analogous fashion to the mediator 
fields in the supersymmetric model given by Eq. ([3]), so this provides a natural framework for warm inflation. During 
inflation, these fields are heavy and overdamped at the origin, yielding a constant cosmological constant that drives 
the accelerated expansion. This flat potential may be lifted by radiative corrections, soft SUSY breaking mass terms, 
non-renormalisable interactions or supergravity corrections, typically yielding a leading correction of the form: 

V = V (i + 1(JlY) , (34) 



n \m p f 

for n > 0, while for n — the dominant contribution comes from logarithmic radiative terms: 

V = Vq ( 1 + 7ln (^)) • (35) 

Here, we will focus on the most common scenarios where either the n = or the n — 2 term dominates, both in the 
strong and weak dissipative regimes. The slow-roll parameters for the hybrid models are given by: 

2 / \ 2 — 2n / \ 2 — n / \ 2 

«* = y(t) • *-<-'>•> ft) • '*-^Ty & = 2(»-D("-2)ft) ,(36, 

where we assumed that, at horizon-crossing, the potential is dominated by the constant term, Vq. This limits the 
values of 7 and 4> and hence the validity of consistency relations between observables. Note that when the corrections 
dominate, the potential behaves effectively as a monomial and the relations derived earlier apply. 

For the n = 2 hybrid model, n s > 1 as for the quadratic monomial potential, in both the weak and the strong 
dissipation regimes at horizon crossing. Since this is disfavoured by the Planck data, we will not consider the 
consistency conditions in this case, focusing only on the logarithmic hybrid potential. 
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1. n = 0; Weak dissipative regime 
For the logarithmic hybrid model in Eq.p5p. we obtain the following consistency relations: 

B B = 13(1 - n s ) ± 8y/2{n s - l) 2 - < , n t = (1 - n s ) ± \/2(n s - l) 2 - < . (37) 

Note that the '+' solutions correspond to 7 < 1, while the ' — ' solutions hold for 7 > 1. We find that only the ' — ' 
solution has a region of parameter space in which the corrections are small and the consistency relation predictions 
are compatible with the Planck results. 

Tensor spectral index 




FIG. 4: Contours of n t in the (n s ,n' s ) plane for the n — hybrid model in the weak dissipative regime, using the ' — ' solution. 
We exclude the grey regions where the consistency relation predicts a positive or imaginary n t . In the blue cross-hatched region 
the radiative corrections dominate, 7I \n((j>/m p )\ > 1. 



Baryon Isocurvature fraction 




0.92 0.93 0.94 0.95 0.96 0.97 0.98 0.99 1.00 



FIG. 5: Contours of Bb in the (n s ,n' g ) plane for the n = hybrid model in the weak dissipative regime, using the ' — ' solution. 
The light grey region indicates where Bb exceeds WMAP bounds. In the blue cross-hatched region the radiative corrections 
dominate, 7|ln(0/m p )| > 1. 

As illustrated in Figures |4] and El the logarithmic corrections dominate in a large fraction of the parameter space, 
where the potential behaves effectively as a single monomial. The constant part only dominates in a narrow region 
of parameter space, where n' s ~ 0.005, n t ~ —0.01 and we find relatively large, positive values of Bb ~ 0.2 — 0.4. 
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2. n = 0; Strong dissipative regime 

For the logarithmic hybrid model, n = 0, in the strong dissipative regime we find the following relation: 

n t = i (2(1 - n s ) ± - l) 2 - 126<) . (38) 

Note that the '+' and ' — ' solutions correspond to 7 < 7/12 and 7 > 7/12, respectively. However, the '+' solution 
only yields nt < for n s > 2 and hence is not consistent within the observationally allowed region. In Figure [5] we 
illustrate the contours of n t for the ' — ' solution. 

Tensor spectral index 

0.04 
0.02 
0.00 
-0.02 



-0.04 



0.92 0.93 0.94 0.95 0.96 0.97 0.98 0.99 1.00 




FIG. 6: Contours of nt in the (n s ,n' s ) plane for the hybrid model with n — in the strong dissipative regime. We exclude the 
grey region where the consistency relation yields nt > or imaginary. 

For this relation, we find that Q 3> 1 is satisfied and the logarithmic corrections are subdominant within the region 
where n t < for n s > 1. For n s < 1 keeping the logarithmic corrections subdominant gives an upper bound on n' s , 
for example for the central value of n s ~ 0.96 we find \n' s \ > 0.005 with n' s < 0. 



C. Hilltop potentials 

Hilltop scenarios are small field models where inflation occurs near a maximum of the potential, being parametrized 
by: 

/ ji/f / j, \ 2 \ 

(39) 

Hence, this is equivalent to the n — 2 hybrid model with a negative mass squared, 7 —¥ —M. The consistency relations 
are thus the same for the hilltop model as for the n = 2 hybrid model, although the requirement M{4>/m p ) 2 /2 < 1 or 
> yields distinct constraints on the parameters. 




1. Weak dissipative regime 



In the weak dissipative regime we derive the following consistency relation for the hilltop potential: 

In' n' 

B B = (l-n s ) '-, nt = — s 40 

n s - 1 2(1 -n s ) 

The only simultaneously allowed regions of parameter space for Bb and nt correspond to n s < 1 with negative 
running, n' s ~ —0.001, as shown in Figure[7] Within this region we find n t ~ —0.01 and ~ —0.05 < Bb ^ 0.1. 
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FIG. 7: Contours of nt (left) and Bb (right) in the (n s ,n' s ) plane for the hilltop model in the weak dissipative regime. The 
consistency relations and exclusions are the same as for the n — 2 hybrid model, although the region in which the calculation 
is valid is distinct. 



2. Strong dissipative regime 



In the strong dissipative regime we derive the following consistency relation: 

nt = ~ (l4(l - n s ) ± ^154(71, - I) 2 - 315n' s ) . (41) 

However, only the ' — ' solution yields regions of parameter space where the constant term dominates and the consis- 
tency relation is compatible with the observational results. 

Tensor spectral index 



"s? 0.000 



/ \ 


■■ P/andr+WP-BAO:ACDM + dnjdhk 

H P/anc*+WP+BAO:ACDM+ dn./dlnk +i 


-v'b ••1 











L 0.92 0.93 0.94 0.95 0.96 0.97 0.98 0.99 1.00 



FIG. 8: Contours of nt in the (n a ,n' s ) plane for the hilltop model in the strong dissipative regime. We exclude the grey region 
where the consistency relation yields n t > 0. The blue cross-hatched region indicates where M((f>/m p ) 2 /2 > 1 or < 0. 

Although the region where the consistency relation is valid is rather small in this case, it is in an observationally 
relevant range. In particular, it predicts nt ~ —0.02 with a small but non-vanishing spectral running, n' s ~ —0.002. 



V. DISCUSSION 



In this work, we have derived consistency relations between observables for warm inflation models based on su- 
persymmetric two-stage interactions in the low-temperature regime. We have focused on the most well-studied case 
where dissipation is dominantly generated by the excitation and decay of heavy virtual modes coupled directly to 
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the inflaton. We have obtained relations amongst the parameters characterizing the spectrum of adiabatic perturba- 
tions in the scalar and tensor sectors, as well as the fraction of baryon isocurvature perturbations and its spectral 
index associated to the generation of the observed baryon asymmetry during warm inflation through the dissipative 
interactions. 

The relations between the different observables are, in general, rather intricate and it is not always possible to 
express them in a simple analytical form. However, we have shown that in the weak and strong dissipative regimes 
one may obtain generic relations that are independent of the inflationary potential and which constitute the main 
result of this work, in Eqs. ([25 ]) . ([25 ]) and ([50 )1 . 

In addition we have noted that the thermal origin of the adiabatic perturbations modifies the relation between the 
tensor-to-scalar ratio and the tensor spectral index, giving r < 8|n t |, which is always smaller than the corresponding 
cold inflation prediction. Although this is also a feature of supercooled models with multiple scalar fields, this 
inequality may be used along with the consistency relations derived in this work to probe the thermal nature of the 
inflationary state of the Universe. 

Prior to this year's release, the bounds on the baryon isocurvature perturbations were expected to increase by 
an order of magnitude and the tensor-to-scalar ratio was expected to be probed down to ~ 0.05 [27]| . Interestingly 
the bounds on isocurvature appear more relaxed compared to WMAP data, but as pointed out in (26| the relaxed 
constraints arise from the observed deficit in the Sachs- Wolfe plateau and not due to the expected phase shift in the 
acoustic oscillations which would be a smoking- gun for an isocurvature component. This effect could have several 
different explanations and is thus not a definite sign of isocurvature modes, although the data allows for a significant 
admixture of adiabatic and isocurvature components. 

In the warm baryogenesis scenario, the correlation between the produced baryon isocurvature perturbations and 
the main adiabatic fluctuations, as well as the isocurvature spectral index, depend on the form of the inflationary 
potential, yielding in general less stringent bounds on the isocurvature fraction than for example in the axion or 
curvaton scenarios. However, when considering specific inflationary potentials the model predicts a spectral index 
that can be used to derive stronger bounds on the isocurvature contribution. The Planck collaboration will complete 
and release its CMB polarization data analysis next year, with a revised analysis of inflationary parameters. One 
then expects the constraints on the relevant observables to increase significantly and hopefully for new observational 
signatures to arise as well. However, even if an isocurvature component is unambiguously found, Planck will be unable 
to distinguish between cold dark matter and baryon isocurvature perturbations, with a promising observational avenue 
being future 21 cm fluctuation surveys (28|. 

The Planck results for the scalar spectral index and associated running are even more constraining. With only the 
WMAP 9-year data the central value of the spectral index is red (blue) without (with) running, although consistent 
with a scale invariant power spectrum in both cases at the 95% C.L. A combination of the WMAP data with eCMB, 
BAO and Hq helps alleviating this degeneracy, at the same time ruling out the Harrison- Zcldovich spectrum at the 
95% C.L. and yielding a negative running 25j. On the other hand, the Planck collaboration has found robust evidence 
for a red spectral index, with and without including running, thus alleviating this degeneracy. This disfavours the 
quadratic monomial and hybrid potentials in warm inflation, given that n s > 1 in both cases. We note that although 
in cold inflation the quartic potential is now ruled out, in warm inflation, on the other hand, it shows an impressive 
agreement with data |5|. 

Warm inflation models may also lead to observable nongaussian signatures in the primordial spectrum from the 
coupling between the inflaton and radiation fluctuations. These depend on the form and strength of the dissipation 
coefficient, as well as on the viscous properties of the radiation bath, although in general we may expect f^L ^ 10 
[29H3lL l3~il [Hf. Planck has so far found no evidence for primordial non-gaussianity, and a dedicated analysis of 
warm inflationary models for a T-independent dissipation coefficient places an upper bound on the dissipative ratio 
at horizon crossing, Q < 400 [26| . which represents an improvement of two orders of magnitude with respect to the 
earlier WMAP bounds. Although this is not directly applicable to the case considered in this work, we expect that 
further exploration of nongaussianity in general regimes of warm inflation combined with the results presented in this 
work may help constraining different models of warm inflation in the future. 

Many of the consistency relations derived in this work require a knowledge of the properties of the tensor power 
spectrum, in particular its (running) scale dependence. In this sense, some of our results may only be tested with the 
next generation of CMB polarization experiments. Nevertheless, we expect Planck's analysis of the CMB polarization 
spectrum to yield much tighter constraints on cold and warm inflation models next year. 

Although our analysis has focused on a particular realization of warm inflation that has been well established in 
the context of quantum field theory, our results should be put in a wider perspective of an alternative paradigm 
where thermal or even non-thermal statistical states describe the evolution of the inflationary universe. In particular, 
we expect other related models to yield relations between CMB observables that can be directly tested against 
observational data, such as those presented in this work. Hence, we hope that our work motivates further research in 
this direction from both the observational and theoretical perspectives. 
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